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Abstract 

We make a first crude direct estimate of the net gluon polarization in the pro- 
ton, AG, combining data on the asymmetries in high-pr hadron production 
from the HERMES, SMC and COMPASS collaborations. Although these 
data sample a restricted range of x, they provide no hint that AG is large. 
Fixing the normalizations of different theoretical parametrizations using the 
hadron asymmetry data, we find typical central values of AG ~ 0.5, with 
uncertainties of similar magnitude. Values of AG > 2 are disfavoured by 
Ax 2 ~ 9 to 20, depending on the parametrization used. 
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1 Introduction 



The spin structure of the proton is still uncertain, despite considerable ex- 
perimental effort, theoretical ingenuity and several surprises jT]. It certainly 
differs from naive expectations formulated within the non-relativistic quark 
model, with strange quarks apparently polarized oppositely to the proton: 
As < 2j, and the quarks altogether apparently contributing only about 
30% of the proton spin: AS = Au + Ad + As ~ 0.3 ±0.1 2J. However, 
whether the remainder of the proton spin is due to gluons AG and/or orbital 
angular momentum L z remains an open question pQ. One particular theo- 
retical interpretation of the proton spin is provided by chiral soliton models, 
which suggest that AS, AG ~ 0, with L z ~ 0.5, up to corrections of higher 
order in 1/N C and m s jSj- 

The question of the magnitude of AG was given a high profile by sug- 
gestions that gluons might be making a significant negative contribution to 
the net quark spins, which might even be sufficient to explain all the 

negative polarization of the strange quarks: As ~ — (a/2n) x AG, which 
would require AG ^2 *. Thus there are two conceptually interesting values 
for AG: a value ~ 1/2 which would be comparable to the proton spin and 
hence make an important direct contribution to it, and a value ~ 2 which 
could 'explain away' the strange quark spin. It is worth keeping in mind that 
both of these effects require positive AG, and that very few models suggest 
a negative value. 

First attempts to estimate AG from data have been made using measure- 
ments of the polarized structure function gi at different momentum scales 
Q 2 to extract AG indirectly using next-to-leading order (NLO) QCD [HIE]- 
These first attempts have been inconclusive, unable to exclude any of the 
interesting theoretical possibilities AG ~ 0, 0.5 or 2. It is hoped that new 
structure function measurements, e.g., by COMPASS jTU], will be able to 
refine this indirect extraction of AG. 

However, the main objective of COMPASS is the direct determination of 
AG via asymmetries in the production of cc and high-pr hadron pairs. The 
direct determination of AG is also a key objective of the RHIC polarized 
beam programme ^T]. First results from both PHENIX at RHIC [12] and 

"However, this suggestion is controversial, in particular because the gluon renormal- 
ization of the quark spin is scheme-dependent . 
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COMPASS at both high jT3j and low Q 2 ^3] have recently been announced, 
the former on the asymmetry All, and the latter on the asymmetry in high- 
er hadron-pair production. The high-Q 2 COMPASS result was in fact the 
third result on this asymmetry, having been preceded by measurements by 
the HERMES [T^ and SMC [T7] collaborations. 

The current PHENIX result is difficult to analyze in terms of AG, since 
the error is still quite large, the central value is outside the physical re- 
gion, and the relation to AG is not single-valued. The errors in the high- 
Q 2 measurement by COMPASS are no smaller than those of the previous 
HERMES and SMC measurements, but the errors in their low-Q 2 are signif- 
icantly smaller, and COMPASS benefits in the comparison with HERMES 
from having more generous kinematics, which reduces potential issues re- 
lated to higher-twist effects. In fact, as we discuss in more detail below, the 
COMPASS, SMC and HERMES data seem to be telling a consistent story. 

In this paper we make a preliminary combination of the HERMES [To] . 
SMC PH and COMPASS PUd data, seeking a first direct indication of the 
possible magnitude of AG based on hadron-asymmetry data. Despite their 
limited precision, these data are already precise enough for such an analysis 
to be carried out within the framework of existing parametrizations of the 
possible polarized gluon distribution. On the other hand, they are not yet 
sufficiently precise to merit a fully-fledged NLO fit. 

We find that AG is unlikely to be as large as was desired in attempts to 
'explain away' the negative value of As, though a substantial gluonic contri- 
bution to the spin of the proton can certainly not be excluded. Using three 
different parametrizations of the polarized gluon distribution, we find central 
values of AG ~ 0.5, with errors of similar magnitude. The suggestion that 
AG > 2 is disfavoured by A% 2 ~ 9 to 20, depending on the parametrization 
adopted, assuming that the polarized gluon distribution does not exhibit 
unexpected behaviour outside the limited x range covered by the current 
experiments. 

2 Available Experimental Information 

We first review the relevant experimental information that is currently avail- 
able. Three recent attempts to extract AG from NLO analyses of deep- 
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inelastic structure function data yield the following estimates: 

AG = 1.026 ± 0.549 Set 3, Ref. = BB3 

AG = 0.931 ± 0.669 Set 4, Ref. = BB4 (1) 

AG = 0.533 ± 1.931 Ref. & = AAC 

The first two estimates were made by the same group, and the difference be- 
tween their central values may be indicative of the systematic errors in this 
indirect approach to AG that are associated with the choice of parametriza- 
tion of the polarized gluon distribution. The larger error in the third estimate 
may mark a more realistic assessment of the systematic errors in this indi- 
rect approach. According to these analyses, each of the parametrizations 
used in jHJ |H] would be compatible with AG = 2 at the 2-cr level. 

The measurements of the high-p^ hadron production asymmetry that we 
use here are the following: 

HERMES US] : AG/G = 0.41 ± 0.18 ± 0.03 0.06 < x G < 0.28, 

SMC [T7j : AG/G = -0.20 ±0.28 ±0.10 (x G ) = 0.07, 

COMPASS HSj : AG/G = 0.06 ± 0.31 ± 0.06 (x G ) = 0.13, 

Q 2 > 1 GeV 2 

COMPASS [H| : AG/G = 0.024 ± 0.089 ± 0.057 (x G ) = 0.095. 

Q 2 < 1 GeV 2 

(2) 

where we have indicated in each case the available information on the kine- 
matic range of the measurement. 

In addition, COMPASS has recently released preliminary results ^S] for 
gluon polarization from open charm, based on the 2002-2003 data: 
AG/G = -1.08 ±0.73 at (x G = 0.15), RMS=0.08. This channel has very lit- 
tle background, but very low statistics, resulting in large errors. We therefore 
do not use these preliminary open charm data in our fits, but note that that 
the central value of AG is negative, providing additional qualitative evidence 
against a large positive AG. 

In order to convert the measurements (j2J) into estimates of AG, one needs 
to assume a suitable form for the unpolarized gluon distribution G(x, Q 2 ) and 
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specify the relevant momentum transfer scale Q 2 . As our defaults, we use a 
recent MRST [THj gluon distribution and assume that Q 2 ~ 5 GeV 2 , and we 
discuss later the sensitivity to the assumed value of Q 2 . 



3 Fits to Asymmetry Data 

In making our direct estimates of Q 2 , we assume three trial forms for the 
polarized gluon distribution AG(x,Q 2 ), proposed by the groups mentioned 
earlier (HI E]. The explicit expressions for the three parametrizations are 
available as FORTRAN codes that can be downloaded from the HEPDATA 
site [19J. Schematically, they can be written in the following form which 
emphasizes the overall normalization: 



We treat the overall normalizations as free parameters, but retain as 
defaults the values of the other parameters chosen in [H] and the central 
fit values found in [§]. We then fit the overall normalizations to the three 
asymmetries (j2J), and hence determine the integrated gluon polarization AG. 

The limited precision of the presently available asymmetry data is insuf- 
ficient to merit a fully-fledged multi-parameter NLO fit. Even with the much 
more precise deep- inelastic scattering data set, the authors of Ref. jH] found 
it necessary to to give up on a fully-fledged fit to the parameters of the form 
they proposed for AG(x, Q 2 ), imposing by hand several constraints and leav- 
ing the overall normalization as the main free parameter. On the other hand, 
Ref. [H] did not impose supplementary constraints on their parametrization 
of the polarized gluon distribution, which is why the errors for AG that they 
quote are considerably larger. 

Our procedure is clearly hostage to unforeseen properties of the polarized 
gluon distribution AG(x, Q 2 ) at values of x outside the experimental ranges 
given in (j2J). This possibility would introduce a systematic error that we are 
unable to quantify. 

Before fitting the values of the normalizations in the different para- 
metrizations, we have first assumed values that correspond to an integral 



BB3 : 
BB4 : 
AAC : 



AG(x,Q 2 ) 
AG(x,Q 2 ) 
AG(x,Q 2 ) 



A B B3-fBB 3 (x,Q 2 )M, 
Abba • fBB&(%, Q 2 ) 0, 
Aaac- fAAc(x,Q 2 )M- 



(3) 
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Figure 1: Experimental results for AG/G [7^1 f"F6| |J7| / compared with the 
theoretical parametrizations 0, using normalizations adjusted to yield 
AG = 2. 



AG = 2 *, and assessed the goodness of fit by evaluating the corresponding 
X 2 functions. As seen in Fig. [TJ each of the parametrizations ® reproduces 
the general trend of the values of AG indicated by the data from the SMC, 
COMPASS and HERMES. However, each of the parametrizations also vio- 
lates the unitarity bound AG(x) < G(x) at larger x. Ignoring this problem 
for the time being, we find x 2 = 16.4, 10.5 and 22.5 for the BB3, BB4 and 
AAC functions, respectively. If one caps each polarized distribution by the 
unitarity bound, keeping the same normalization at lower x where there are 
measurements, the integrals are reduced to AG = 1.63, 1.99 and 1.80, re- 
spectively. If one were to attempt to compensate for these reductions in AG 
by increasing the normalization factors Abb3,bb4,aac, t ne x 2 values would 
each be increased. 

Much better fits to the parametrizations (j3J) are obtained when the over- 
all normalizations are allowed to float: x 2 = 1-1, 1-5 and 2.5 and, as seen in 
Fig. El each parametrization now respects the unitarity bound AG(x) < G(x) 

*We repeat that such a value was compatible with the NLO analyses of jEHU], though 
not the central value suggested by their fits. 
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Figure 2: Experimental results for AG /G fTh^. \l% f73f [Tffi compared with the 
theoretical parametrizations ' L 8, 9], with normalizations adjusted to yield 
best fits to the data. 



for all x. In view of all the assumptions and uncertainties, we are reluctant 
to quote the differences in % 2 , A% 2 = 15.3, 9.0 and 20.0 , as numbers of 
standard deviations by which AG > 2 is disfavoured. In particular, the A% 2 
found for the AAC fit would change significantly if the full freedom of the 
parametrization were explored. However, it is clear that the current asym- 
metry data offer no hint in favour of the option that AG > 2. The following 
are the best fit values of AG that we find for each of the parametrizations 
(J3J) and the formal errors: 

BB3 : AG = 0.31 ±0.43, x* = 1-1 , (4) 
BB4 : AG = 0.39 ±0.54, X 2 = 1-5 , 
AAC: AG = 0.57 ±0.32 X 2 = 2.5. 

The best fits are compared with the data in Fig. 121 

In the preceding discussion we have compared the experimental data with 
polarized gluon parametrizations AG(x, Q 2 ) evaluated at a canonical value 
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Q 2 = 5 GeV 2 . However, each of the measured values of AG/G in eq. (J2J) 
results from averaging by the relevant experiments over a range of values 
of x and Q 2 . The limited information in the experimental papers and the 
low precision of the currently available data make it impossible to provide 
an accurate estimate of systematic errors due to this averaging. We can, 
however, get a semi-quantitative estimate of the relevant error by repeating 
the fitting procedure at several values of Q 2 , as shown in Table I. 



Table I 





BB4 


BB3 


AAC 


Q 2 




x 2 


x 2 




x 2 


x 2 




x 2 


x 2 


GeV 2 


AG 


best 


for 


AG 


best 


for 


AG 


best 


for 






fit 


AG=2 




fit 


AG=2 




fit 


AG=2 


1.5 


0.31±0.40 


1.6 


19.7 


0.46±0.44 


1.0 


13.4 


0.58±0.31 


2.2 


22.8 


2.0 


0.33±0.43 


1.5 


16.5 


0.33±0.37 


1.2 


21.2 


0.57±0.31 


2.3 


23.4 


5.0 


0.39±0.54 


1.5 


10.5 


0.31±0.43 


1.1 


16.4 


0.57±0.32 


2.5 


22.5 


10.0 


0.43±0.62 


1.5 


6.5 


0.30±0.48 


1.1 


13.5 


0.58±0.33 


2.6 


20.6 



Fits to the HERMES, SMC and COMPASS AG/G data 0J, forQ 2 = 1.5, 2, 5 and 
10 GeV 2 , using the parametrizations 0] of G(x,Q 2 ). For each parametrization 
we list the best- fit value of AG and its x 2 > as we H as the x 2 value corresponding 
to AG = 2. 

The trend shown by these fits is clear: for 1.5 < Q 2 < 10 GeV 2 (which 
includes the preferred value Q 2 = 3 GeV 2 quoted in [14J and the three 
parametrizations J3J), the best fit values AG range between 0.30 ± 0.48 and 
0.58 ± 0.33. In all cases the AG = 2 value is significantly disfavoured. 

The current discussion is based on the COMPASS 2002 - 2003 data with 
Q 2 > 1 GeV 2 ^3] and approximately 10 times more data with Q 2 < 1 GeV 2 
|14j . The present statistics will be approximately doubled with the 2004 data, 
When the full 2004 set of all Q 2 COMPASS data is analyzed, the statistical 
error on AG/G is expected to go down to ±0.05 [TSJ, compared with ±0.31 
in the currently available 2002-03 data set. However, the inclusion of Q 2 < 1 
GeV 2 data introduces additional theoretical uncertainties. 
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Still, it is clear that the inclusion of the full 2004 data set will significantly 
increase the precision with which the first moment AG can be estimated. To 
see the effect of the increased precision of the future COMPASS data, we 
have repeated the current fits, setting the COMPASS AG/G error at ±0.06, 
while keeping the HERMES and SMC data unchanged. The expected error 
on AG shrinks down to ±0.44, ±0.55 and ±0.16, for the BB3, BB4 and AAC 
parametrizations, respectively. In the cases of the two BB parametrizations, 
the expected error reduction is rather modest, but we attribute this to the 
fact that in this exercise the fit includes only one x value with increased 
precision. This underlines the importance of providing the high-precision 
values of AG/G over a range of x. 



4 Summary and Prospects 

We have made in this paper a first direct estimate of the net gluon polar- 
ization in the nucleon, based on hadron-asymmetry data in deep-inelastic 
scattering [T3J EH El EI] ■ Despite being very crude and incomplete in its 
kinematic coverage, this direct estimate has an error that is comparable with 
that provided indirectly by NLO analyses of deep-inelastic structure func- 
tions [HI Ej. We find a favoured value of AG ~ 0.5, with a formal error of 
similar magnitude. Values of AG > 2 are disfavoured by A% 2 ~ 9 to 20, 
depending on the parametrization of the polarized-gluon distributions that 
is used. 

There are good prospects for a significant improvement soon in the ac- 
curacy with which AG is known, thanks to new data from COMPASS and 
RHIC. The present data are insufficient to exclude strongly the hypothesis 
that all the apparent negative value of As might be induced by gluons via 
renormalization in one particular scheme. The forthcoming data should be 
able to resolve this issue. However, they might not be able to determine 
whether gluons carry a large part of the nucleon spin, AG ~ 1/2, or whether 
their contribution is as small as that due to the quarks, as expected quali- 
tatively in chiral soliton models j^j. There are surely still many interesting 
twists and turns still to come in our understanding of the nucleon spin, but 
direct determinations of the gluon spin now seem poised to make an impor- 
tant step forward. 
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